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Abstract 

> 

O. 

Using both the Born-Oppenheimer idea and the de Broglie-Bohm interpreta- 

O 

tion of wavefunction we represent in a different way the semiclassical quantum 
, gravity from the Wheeler-DeWitt equation in an oscillating regime which can 

q-i preserve completely the unitary quantum evolution of a matter field at the 

(3JT), expense of a nonlinear gravitational field equation, but has the same asymp- 



X 



totic limit as the others. We apply the de Broglie-Bohm interpretation to the 



nonlinear gravitational field equation to develop a perturbation method to 
find the quantum corrections of a matter field to the gravity. The semiclas- 
sical Einstein equation with the quantum corrections is found for a minimal 
quantum FRW cosmological model. 
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I. INTRODUCTION 



The problem of unitarity and back-reaction of a matter field in a curved spacetime has 
been an important and critical issue of gravity such as black holes and the very early Uni- 
verse. Although the complete solution of this problem should be found within the context of 
the right and consistent quantum gravity, there has been quite recently an attempt to tackle 
this problem from the canonical quantum gravity based on the Wheeler-DeWitt (WDW) 
equation ffl. In the canonical quantum gravity approach to the gravity plus matter sys- 
tem, the gravitational wavefunctions of the WDW equation are peaked along (semi) classical 
trajectories of the (semi) classical Einstein equation, and the matter field, which obeys a 
time-dependent functional Schrodinger equation, evolves quantum mechanically along each 
trajectory These two equations constitute the so-called semiclassical quantum gravity. We 
shall distinguish this approach from the conventional field theoretic approach in the curved 
spacetime 0] . Recently the semiclassical quantum gravity has been applied to a two dimen- 
sional black hole model [||]. 

One can derive the semiclassical quantum gravity in many different ways, for instance, 
depending on whether or not he or she includes the back-reaction of the matter field to the 
semiclassical Einstein equation in the lowest order of the power series expansion in terms of 
the inverse Planck mass. When one searches for a WKB-type wavefunction for the entire 
WDW equation and expands its total action in the inverse power series of the Planck mass, 
he or she obtains the vacuum (classical) Einstein equation and the time-dependent functional 
Schrodinger equation for the matter field on it Q. On the other hand, when one adopts 
the Born-Oppenheimer idea by treating the gravitational field as a massive particle degree 
of freedom and the matter field as a light particle degree of freedom, he or she can separate 
two different mass scales and derive the semiclassical Einstein equation with the quantum 
expectation value of the energy-momentum tensor operator as a source of matter and the 
time-dependent Schrodinger equation for the matter field [0,0. In previous papers 0, by 
separating and expanding the WDW equation in the inverse power series of the Planck mass, 
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we were able to derive the semiclassical quantum gravity and to find the quantum states of 
the matter field in the asymptotic limit of large Planck mass. Quite recently Bertoni et al 
H have shown that the semiclassical quantum gravity represented in three different ways 
are in fact equivalent to each other in the asymptotic limit Oi^^j. 

In the asymptotic limit of large Planck mass, the semiclassical quantum gravity is uni- 
tary in the sense of conservation of probability of the matter wavefunctions 0. In a dif- 
ferent representation, there are, however, unitarity violating terms in the time- dependent 
Schrodinger equation for the matter field beyond the asymptotic limit |§. These unitarity 
violating terms come from the slowly varying amplitude of the gravitational wavefunction 
in the expansion of the WDW equation in the inverse power series of the Planck mass. The 
origin of these unitarity violating terms may not, however, follow inherently from the WDW 
equation. Bertoni et al || have further shown that the semiclassical quantum gravity can 
indeed preserve the unitarity of quantum fields provided that one can define the cosmolog- 
ical time defined through the Hamilton- Jacobi equation of gravitational field including the 
back-reaction of matter field. A remark on the different representations of almost the same 
semiclassical quantum gravity is that the semiclassical quantum gravity in Ref. ]7[ is based 
on the complicated matrix equation of the WDW equation which incorporates actually the 
scattering of Cauchy data at the second or third quantized level and also considers linear 
combinations of the wave functions | lQfl , in contrast with one wave function in Refs. |J and 



||. So, it would be worthy to show the unitarity of quantum fields from the WDW equation 
using the equivalent matrix equation 

In this paper we observe that the WDW equation is not uniquely separated into the 
gravitational field equation and the time-dependent Schrodinger equation according to the 
different mass scales following only the Born-Oppenheimer idea 0, and that, especially when 
we apply de Broglie-Bohm interpretation ]TT] at the same time to the wavefunction of the 
gravitational field equation, this arbitrariness in the separation of the WDW equation in an 
oscillating regime can also allow the semiclassical quantum gravity which consists of a non- 
linear gravitational field equation and a time-dependent Schrodinger equation that indeed 



does preserve the norm of the quantum states of matter field during the quantum evolution. 
As in Ref. ||, we assume no asymptotic limit to derive the semiclassical quantum gravity. 
However, in the asymptotic limit, it reduces identically to the semiclassical quantum gravity 
in previous papers 0,0. Thus, the question of unitarity of the quantum evolution of matter 
fields may not be settled within the canonical quantum gravity, and should be answered by 
more fundamental theory. Furthermore, we shall use this semiclassical quantum gravity to 
find the higher order quantum corrections of the matter field to the gravity. Finally, we 
shall apply the semiclassical quantum gravity to the quantum Friedmann-Robertson- Walker 
(FRW) cosmological model coupled to a minimal scalar field. 

The organization of the paper is as follows: In Sect. II we derive the nonlinear gravi- 
tational field equation and the time-dependent Schrodinger equation, and show that they 
preserve the unitarity of the quantum evolution of matter field. In Sect. Ill we develop a 
perturbation method for the semiclassical Einstein equation with the quantum potential and 
the higher order quantum corrections of matter field. In Sect. IV we apply the semiclassi- 
cal quantum gravity to the quantum FRW cosmological model to find the effective energy 
density. 



II. NONLINEAR GRAVITATIONAL FIELD EQUATION AND UNITARITY 



We shall follow the Born-Oppenheimer idea to separate the WDW equation into the 
nonlinear gravitational field equation and the time- dependent Schrodinger equation. The 
canonical quantum gravity for a gravity coupled to a matter field, typically represented by 
a scalar field, is described by the WDW equation 



h 1 



-V 2 - m P V{h a ) + H 



i 5 



4>, h a 



*{h a ,<f>) = 0, 



(1) 



2mp ' ' \hi 

where h a represents the gravitational field and <fi the matter field. Here V 2 is the 



Laplace-Beltrami operator in superspace with the signature ( — ,+, 
Oppenheimer-type wavefunction for the WDW equation 



+). Take the Born- 
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\*(h a ,cf)))=i>(ha)m4>,h a )), (2) 

where ip depends only on the gravitational field and $ both on the matter field and on the 
gravitational field as an explicit parameter. The Born-Oppenheimer idea assumes that the 
quantum states of matter field belong to a Hilbert space that varies continuously on the 
superspace of the gravitational field. The validity of the form of wavefunction (^) should be 
checked by the unitarity of the quantum state of matter field in the semiclassical quantum 
gravity. 

Before we delve into representing in a different way the semiclassical quantum gravity 
that preserves the unitarity of the quantum matter field, we point out a major similarity and 
a minor difference among the approaches to the semiclassical quantum gravity in Refs. ||[7] 
and ||. The similarity is that, as explained in the Introduction, most of the approaches [0-|8| 
take into account the expectation value of energy-momentum tensor as its matter source in 
the semiclassical Einstein equation. This is the quintessence of the Born-Oppenheimer idea 
originally applied to a molecular system, but in Ref. the back-reaction was excluded at 
the lowest order by expanding the total action in the inverse power series of the Planck 
mass. 

The difference comes from the interpretation of the WDW equation (p]) whether as the 
zero-energy Schrodinger equation for a coupled system or as a functional wave equation 
describing the evolution (scattering) of Cauchy data from an initial spacelike hypersurface, 
usually assumed to be the early stage of the Universe, to a later spacelike hypersurface 
[10|1 • 111 a toy model without a matter field, the WDW equation may apparently look like 



the Schrodinger equation with zero energy, but with additional inhomogeneous gravitational 
degrees of freedom or matter field degrees of freedom, Eq. ([!]) has the Minkowski signature 
(— .+, ■ • ■ , +). The exact wave functions for the quantum FRW cosmological model mini- 
mally coupled to a scalar field with a general potential, which is to be considered in Sec. IV, 
are extremely difficult to find due to the non-separability of quantum states of the scalar 
field from the gravitational field and have not been found yet. In a systematic way we can 
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make use of the initial value problem for Eq. (JJ) formulated in a matrix form to describe the 
evolution (scattering) of Cauchy data from one spacelike hypersurface to another, including 
the evolution of any quantum state. The matrix gravitational field equations in Ref. |J in- 
deed incorporated the Cauchy data. The semiclassical quantum gravity in previous papers 
|7]] was based on this exact matrix formulation for the canonical quantum gravity, and took 
into account not only single wave function but also a linear combination of wave functions. 

Following Refs. we may expand any quantum state with respect to a suitably chosen 
orthonormal basis relevant to the matter field: 

=£c fe (/l )|$ fe (Ma)>. (3) 

k 

Substitute Eqs. (Q) and (|3]) into the WDW equation (JJ) to get the matrix equation [0] 

h 2 

a 

( ti 2 

A n /c ( h a ) c k ( h a 



+ ^( h a)(j2 Hnk(h a )Ck(h a ) - - ^nk(h a )c k (h a fj j = 0. 

k^n P k / 



(4) 



where 



Hnk(h a ) : = H|$ fe ) , 
A nk (h a ) : =i($ n |V|$ fc ), 

^nk(K) ■ = V 2 5 nfe - 2iA nk • V + ($ n | V 2 \$ k ) . (5) 

This is an intermediate step for the derivation of the semiclassical quantum gravity. 

We now derive the semiclassical quantum gravity in a different way from the previous 
papers 0, so that it can preserve the unitarity of quantum states completely. A complex 
solution of the gravitational field equation in an oscillating regime can always be rewritten 
in the form 

ij(h a ) = (r^exp^S(h a )y (6) 



Substituting the complex solution (|6|) into the matrix equation (||) and assuming c n ^ 0, we 
obtain the two coupled equations: 



h 2 

-V 2 - mpV(h a ) + n nn (h a ) 



2mp 

h 2 fVc n (h a ) „ A ,,.Ck{h a ) 



/( -rn-'( — 7TT- ~ A nk{h a Hrr ) )i> (ha ) = (7) 



and 



i — VS(/i„) • Vcn(h a ) + — VS(h a ) ■ A nk (h a )c k (h a ) 
mp mp Y 

ft 2 

- V H nfc (/i a )c fc (/i a ) + Qnk{ha)c k {h a ) = 0. (8) 

k+n 2m P k 

Eq. (^D may be interpreted as the nonlinear gravitational field equation with the quantum 
back-reaction of the matter field. In terms of the cosmological time 

JL := — VS(h a ) ■ V, (9) 
or mp 

Eq. (|§D can be rewritten as 

ih—c n (h a ) = V B. nk (h a )c k (h a ) 
or f~f 

h h 2 

VS{h a ) ■ V A nk (h a )c k {h a ) - Vr2 nfc (/i a )c fc (/i a ). (10) 

m P V 2m ^ * 

Eq. (pL0|) may now be interpreted as the time-dependent functional Schrodinger equation for 

the matter field on the spacetime determined by Eq. (]7|) . It should be noted that there can 

be some arbitrariness in separating the WDW equation. Instead of H nn in Eq. (0) one may 

use a quantum back- reaction K. nn (h a ). Then, the nonlinear gravitational field equation (^) 

is replaced with K nn . We rewrite the matter field equation by a vector notation 

ih^-c(h a (r)) = M(h a (r))c{h a )(r)) (11) 
or 

where M is the hermitian matrix 

M(h a (T)) = M (h a (r)) + M 1 (h a (r)), (12) 
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consisted of 

M (h a (r)) = J2 Hnk(h a (r)) - —VS(h a (r)) ■ £ A nk (h a (r)) - ^— ]T Q nk (h a (r)), 

k^n m P k lm P k 

yix(h a (r)) = K nn (h a (r)) - H nn (/i a (r)). (13) 

The unitarity of the quantum matter field is always fully respected 

c t( r ) • c(r) = 1. (14) 

A minor difference of this paper from other related works [£|,[5],|| is that we followed 
closely Ref. [10], which interpreted the WDW equation as the quantum scattering of the 
Cauchy data including the quantum states of matter fields. The difference from previous 
works |7j lies on the fact that we made use of the arbitrariness in order to separate the 
nonlinear gravitational field equation (^) from the WDW equation and the time-dependent 



Schrodinger equation (10) that preserves the unitarity of the quantum evolution of matter 
field. The unitarity of quantum field in the oscillating regime has been first proved by 
Bertoni et al M using the semiclassical quantum gravity in Ref. []5J . 

III. QUANTUM CORRECTIONS 

We apply the de Broglie-Bohm interpretation of wavefunctions to the nonlinear gravi- 
tational field equation (0). The main idea of the de Broglie-Bohm interpretation is that a 
time-dependent Schrodinger equation is mathematically equivalent to a pair of equations, 
the time-dependent Hamilton- Jacobi equation with an additional quantum potential and 



the equation for the conservation of probability [11]. The concept of trajectories with the 
quantum potential provides us with all the same physical predictions as the standard quan- 
tum mechanics. Thus, the de Broglie-Bohm interpretation is another way of interpreting 
quantum mechanics, not a kind of approximation methods such as the WKB method, even 
though the de Broglie-Bohm interpretation has the WKB as an asymptotic limit, when the 
quantum potential is negligible. 



S 



According to the de Broglie-Bohm interpretation, we separate the real and the imaginary 
parts of the nonlinear gravitational field equation. Thus, by substituting / = (il)*%lj) l l 2 and 
equating the real and the imaginary parts of Eq. (|7|), we obtain 

1 {VSf - m P V + H wn - f V l F - — Re(Q nw ) = 0, (15) 



2m p 2m p F mp 



I V F 

-V 2 S + — -VS + Im(Q nn ) = 0, (16) 



where 



VF /Vc n c fc x 
Qnn = -^r-{- i}_^A nk —). (17) 

' <T( k ''H 

If there is no parametric-coupling of gravity with matter, the following term 

V - ^ V ' F (18) 
V quant - - 2mp F (18) 



is the quantum potential, and Eq. ( ]T6|) is the equation for the probability conservation [] I 
The last terms in Eqs. (|T5|) and ([IB] ) are the quantum back-reaction of matter field to the 
gravity. 

In the asymptotic limit O^^-V the matter field equation becomes 

m^\h a ) = ]T n nk (h a ) c f\h a ) - -^vs(°\h a ) ■Y.^nk{h a )cf{h a ) (19) 

and the gravitational field equation reduces to the Hamilton- Jacobi equation of the zeroth 
order 

— (VS^Y - m P V + H„ n = 0. (20) 



mp 

As an orthonormal basis, we choose the exact quantum states of the time-dependent func- 
tional Schrodinger equation 

^A|$(o) ( ^ r)) = H(lA ^h a {r))\^\^r)). (21) 

This orthonormal basis simplifies considerably the algebra; especially, the matter field equa- 
tion becomes 



ihf$\h a ) = —LvsV>(ha) ■ A nn (h a )$\h a ). (22) 
or mp 

The solution is 

ci 0) =exp{iJ(A nn ) a dh a ). (23) 

To find the gravitational field wavefunction and the quantum states of matter field, we 
use a perturbation method that expands the coefficient functions and the gravitational action 

as 

oo * 

Cn(/0=£(-f)\ W (U 
p =0 m P 
oo * 

c k {ha) = Y,{— ) P #\ha), (k^n), (24) 

p=l m p / 



and 



~ x Ulp ' 



=0 



oo t 

nha) = E{-^Tf {p) (h a ). (25) 

The dominant contribution to the quantum state comes from c n . The other terms (A; ^ n) 
are of the order of O(^). The prefactor / of the gravitational wavefunction is determined 



by 



whose solution is 



^•V^) = -iv 2 ^°), (26) 



f(o) = \ (27) 

(V5(°)-V5(°)) 1/4 ' 



The first order coefficient functions of the quantum state are 



i j-c£\h a ) = - 1 -{l kn (ha)c i °\h a ), (28) 



where 
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^ kn{ha) {h a ) — ^2A n fc • A nri A nri • Afc n (A • A) nn Jc^(/i ). 



(29) 



We determine S* 1 - 1 ^ from 




2 /(°) 



(30) 



The procedure can be repeated to yield the higher order quantum corrections of the matter 
field to the semiclassical Einstein equation and the exact quantum state. 



We now apply the semiclassical quantum gravity developed in this paper to the frequently 
employed quantum FRW cosmological model. The higher order quantum corrections in 
Sect. Ill to the asymptotic semiclassical quantum gravity have not yet been found explicitly 
even for the FRW model, although higher order quantum corrections were found in the 
semiclassical quantum gravity whose lowest order equation is the vacuum Einstein equation 
and which does violate the unitarity at the first order ||. Moreover, since the quantum 
field preserves the unitarity throughout the evolution, the application to the FRW model is 
expected to be particularly useful in searching for the quantum effects of matter field to the 
gravity. 



The simplest but nontrivial FRW quantum cosmological model is with a minimal scalar 
field (inflaton). The WDW equation is 



IV. QUANTUM FRW COSMOLOGICAL MODEL 



A. Minimal Scalar Field 



27rft 2 d 2 
3mpa da 2 



3mp 



V(a) + nQ— ,0,o) *(a,0)=O, 



(31) 



where a is the size of the universe, and 



V(a) = ka- Aa 3 



(32) 
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is the gravitational potential consisting of the three-curvature and the cosmo logical constant, 
and <p denotes the massive scalar field. The extended supermetric is 



ds 2 



-ada + a 



3 ji2 



(33) 



and the rescaling a = ^ a recovers the WDW equation of the form in Eq. ({]]). Remem- 
bering the superspace signature (— ), the cosmological time is related to the gravitational 
action by 

4tt dS(a) d 



d_ 



3rripa da da 



Then, we find that 



d(r) 



4-7T 



-S' 



(34) 



(35) 



3mpa 

where and hereafter overdots and primes will denote derivative with respect to r and a, 
respectively. We rewrite the nonlinear gravitational field equation as 



27TH 2 d 2 



3m i 



■V(a) + H r 



4vrr F' ,C 



3m P ada 2 8ir v ' "" 3mpa F ^c n ' An S' *t" <' 



ip(a) = 0. (36) 



In terms of the cosmological time, the time-dependent functional Schrodinger equation is 
rewritten as 



or 



H„ fc (r)c fe (r) - ft^B„ fc (r)c fc (r) 

k^n k 

2nfi 2 x , ^ 

}_^n nk (T)c k {T), 

3m P a ^ 



(37) 



where 



H nfc (r) 
B nfc (r) 

^nfc(r) 



($„|H|$ fe >, 

i($ n \-ir\Qk), 

or 



a 1 



2 



aar / or 



a 9 



(38) 



- ar^ aor / ar a ar 

In the above equations we made use of the relations ([34]) and (|35|) . Following the de Broglie- 
Bohm interpretation, we substitute the wavefunction ip = Fe% s and separate the real and 
the imaginary parts of the gravitational field equation to obtain 
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2tt „ /2 , 3mp T ^ N „ 4tt^ F 2^ 2 F" 



and 



-5"' + -Via) = H n „ - 

3rripa 871 3mpaa F 



F' 1 TiF' 



3rripa F 



(39) 



(40) 



where 



R r 



(41) 



and Q nn = -p R„„, We solve Eq. 



and let 



F 
"F 



\aa) 



2 ad2 + ^Im(R, 



(42) 



where we again used (|4]) and (pop. Using again ( |35| ) and (f42|), we rewrite Eq. (|39| ) as 



/d\2 1 
("J +- ya 



8tt 



3mpa 3 



H r 



47r/f 



. U Jiri R.e(R nri ) -|- 



1 • 



a 



(43) 



Smpaa Smpa 

We can interpret Eq. ( f23"D as the semiclassical Einstein equation with the quantum back- 
reaction of the matter field and the higher order quantum corrections from the fluctuations 
of matter field and geometry itself. In this sense, 



. U nn R.e(R nri ) 4 



2nh 2 



U + -U 

a 



(44) 



a . — 10 ft "\ till I ' o 1 

Snipaa Smpa 

is the effective energy density from the quantum fluctuation of matter fields. Eq. ( j3~7l) is 
the time-dependent functional Schrodinger equation for the matter field on the spacetime 
determined by Eq. ([13]) . 



B. Massive Scalar Field 

We work out explicitly the specific case of a massive scalar field. The Hamiltonian of the 
massive scalar field is 
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In Ref. 0] a Fock space-representation of quantum states satisfying the time-dependent 
functional Schrddinger equation was constructed in terms of the annihilation and creation 
operators 



S(r) = #(t)tt, - a\r)4>:{r)l b\r) = (6(r)) t 



(46) 



where C is a complex solution of the classical field equation 

Mr) + 3^f\<Pc{r) + m 2 4> c {T) = 0. 
a(T) 

The effective gauge potential is found in an operator form: 



(47) 



B = a(Tp{r)b(T) + (3(r)b 2 (r) + P*{t)P 2 {t), 



(48) 



where 



a(r) = ^a 3 (r)(0*(r)0(r) +m 2 0*(7-)0(r)) 
P(t) = -— I 



(r) + mV(r)) 



(49) 



Since c„ is of the order 0(1) but Cfc(A; 7^ n) of the order O(^)) we obtain the coefficient 
function 



ci 0) (r) =exp(z | B nn (r)* 



The lowest contributions to R nn and XJ nn are found 



Rl°2 = 



nn "> 



U(o) 



1 (ad) 



2 ad 2 ' 

Thus, the semiclassical Einstein equation at this order becomes 

27rh 2 



(-) 2 +Vw 



57T 



3mpa 3 



3mpa 



nn 1 x 



where 



(50) 



(51) 



(52) 



H nn = fta 3 (r)(W -)(0*(r)0(r) + mV(r)0(r)) 



(53) 
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is the quantum expectation value of the Hamiltonian operator. It should be remarked that at 
the lowest order the quantum effects enter the semiclassical Einstein equation only through 
the quantum potential and there is no contribution from the effective gauge potential. The 
contribution from the quantum potential is still of the order 0\-^j, so we recover the 
asymptotic semiclassical Einstein equation in previous papers [J/J. 
The first order correction can be found from the equations: 



s L nn<- n ~ • o 



,(o) 



(54) 



and 



s L nk^ n 



2B„i.B 



/ • a 

(B n fc + -D n k 

V ft. 



,(o) 



(55) 



Integrating the first order matter field equation 



. d 



2%h 
3a 



ft, r<°> 



(56) 



we get 



c(o) + -MW 

(!) u n ^ my'" 



R_ r (l) 



,(o) 



R J 1 ) 



nfc" 



'> + rn P C k 



_s_ r (i) 

mp 



(57) 



By substituting Eq. (|57|) into Eq. (|43|), we obtain the semiclassical Einstein equation up to 



the first order in — . 

mp 



V. CONCLUSION 

The unitarity of quantum field in a curved spacetime is an important issue in quan- 
tum gravity. In this paper, by applying both the Born-Oppenheimer idea and the de 
Broglie-Bohm interpretation, we derived in a different way the semiclassical quantum grav- 
ity which consists of the nonlinear gravitational equation and the time-dependent functional 
Schrodinger equation that does preserve the unitarity of quantum field without assuming any 
limit. In particular, we applied the de Broglie-Bohm interpretation to the nonlinear gravi- 
tational equation, whose real part is nothing but the Hamiltonian- Jacobi equation with the 
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quantum potential and the contributions from the effective gauge potential of matter fields. 
In an oscillating regime of gravity, we were able to develop a perturbation method for the 
semiclassical Einstein equation, which include the quantum corrections of the matter field. 
Finally, we applied the perturbation method to the quantum FRW cosmological model with 
a minimal scalar field (inflaton) and obtained the semiclassical Einstein equation P3|). The 
right hand side of the equation, the effective energy density, is the sum of the expectation 
value of the Hamiltonian operator, the higher order quantum corrections from the matter 
field and the quantum potential coupled to the effective gauge potential. Since the quantum 
corrections have the order of 0(^-^j, the semiclassical Einstein equation reduces to the 
same asymptotic semiclassical Einstein equation obtained already in Ref. J7|. 

This paper does not, however, resolve completely the problem of unitarity of quantum 
field in quantum gravity and leave it to a more fundamental theory. As pointed out earlier, 
bearing the de Broglie-Bohm interpretation in mind, we made use of the freedom in the 
separation of the Wheeler-DeWitt equation into the gravitational field equation and the time- 
dependent Schrodinger equation, and derived explicitly the semiclassical quantum gravity 
that does indeed preserve the unitarity of quantum evolution of matter field at the expense 
of the nonlinearity of the gravitational field equation. 
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